Noncollinear magnetic order in Quasicrystals 
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Based on Monte- Carlo simulations, the stable magnetization configurations of an antiferromagnet 
on a quasiperiodic tiling are derived theoretically. The exchange coupling is assumed to decrease ex- 
ponentially with the distance between magnetic moments. It is demonstrated that the superposition 
of geometric frustration with the quasiperiodic ordering leads to a three-dimensional noncollinear 
antiferromagnetic spin structure. The structure can be divided into several ordered interpenetrat- 
ing magnetic supertilings of different energy and characteristic wave vector. The number and the 
symmetry of subtilings depend on the quasiperiodic ordering of atoms. 
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The last few years have shown a boom in investigations 
of the spin order in antiferromagnetic films [l|, |3 moti- 
vated by the dramatic changes in the magnetic proper- 
ties of such systems induced by frustration. In contrast 
to the rather well studied spin structure of antiferromag- 
nets on periodic lattices, the antiferromagnetic ordering 
of quasicrystals is subject of ongoing scientific debate. 
Whereas an experimental finding of long-range antiferro- 
magnetic order in rare-earth icosahedral quasicrystals |^ 
turned out to be an artefact 4], theoretical models that 
deal with magnetism on quasicrystals 5] are known to 
exhibit long-range magnetic order. Recent inelastic neu- 
tron scattering experiments on the Zn-Mg-Ho icosahedral 
quasicrystal 6] revealed a very peculiar diffuse scattering 
pattern with icosahedral symmetry at temperatures be- 
low 6K. Such a pattern, in principle, can originate from a 
noncollinear spin arrangement first suggested by Lifshitz 
from pure geometrical considerations P^BS- However, 
real-space magnetic configurations leading to those long 
wave vector correlations remain obscure despite recent in- 
teresting results for quantum spins ^5]. Thus, the knowl- 
edge about the spin structure on quasiperiodic tilings is of 
basic importance for experiments as well as for theoreti- 
cal predictions of new phenomena, which can be expected 
due to nontrivial frustration effects |l3 • 

The patterns found in our theoretical study provide an 
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FIG. 1: Configurations for a frustrated Ising antiferromag- 
net on (a) elementary tiles and (b) six local environments of 
the Ammann-Beenker tiling. Bold lines denote the frustrated 
bonds. The open and filled circles represent different spins. 



explanation for the origin of the antiferromagnetic mod- 
ulations observed experimentally in Ref. 6]. While the 
spin order in antiferromagnets is usually characterized 
by a periodic modulation described by wave vectors on 
the order of inverse atomic distances, the spin order in 
antiferromagnetic quasicrystals admits three-dimensional 
noncollinear structures consisting of several interpene- 
trating subtilings with longer wave vectors. Here we 
report on the details of the low-temperature antiferro- 
magnetic ordering and the map of the local frustration 
for the octagonal tiling. 

We discuss the antiferromagnetic Hamiltonian 
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where S^ is a three- or two-dimensional unit vector in the 
case of classical vector or x?/-spins, and Sf is equal to ±1 
in the case of Ising spins (so Sf = S^ = 0); (z, j) denotes 
the nearest neighbor pairs. For an antiferromagnetic sys- 
tem, the exchange parameter Jij is positive, and neigh- 
boring antiparallel spins contribute a lower energy than 
parallel neighbors. The coefficient Ki is the first-order 
anisotropy constant. Our Monte-Carlo simulations have 
been carried out on finite Ammann-Beenker tilings with 
free boundary conditions. The procedure is a simulated 
annea ling method with at least 15 successive temperature 
steps [l^- At each temperature, the convergence of the 
relaxation process towards equilibrium has been observed 
for any initial configuration after a few thousand Monte 
Carlo steps per spin. Hence, the single-spin-update algo- 
rithm is efficient in our case. At the end of the cooling 
down process, the total energy is just fluctuating around 
its mean equilibrium value. To reduce boundary effects 
only the core of a tiling has been analyzed. The samples 
on the octagonal Ammann-Beenker structure, which we 
shall concentrate on in what follows, are circular, con- 
taining 2193, 11664 and 53018 magnetic moments. 

The octagonal tiling consists of two motifs: a square 
and a rhombus of equal edge lengths a (Fig. ^a)). The 
diagonal bonds are, usually, neglected in the calculations 
[3, [III . We find this disregard physically questionable as 
the exchange coupling increases exponentially with de- 




FIG. 2: The frequency distribution of the energy per spin on 
the octagonal tihng for (a) Ising and (b) vector spins. Sohd 
Hues correspond to the case Jd < 2 J, dashed hues to Jd > 2 J. 
Purely antiferromagnetic interaction at kT — 0.0 IJ is consid- 
ered. Top- views of portions of Monte- Carlo configurations 
with underlying tilings are shown as insets. The light and 
dark circles represent different spins in (a) and different ener- 
gies in (b), respectively. 



creasing interatomic distance. In the present investiga- 
tion, the short diagonal of the rhombus and the sides 
of the motifs have been considered as nearest neighbors. 
We distinguish the two cases Jd > 2J and Jd < 2 J, 
where Jd denotes the interaction along the short diag- 
onal and the interaction strength along the sides J is 
unity. The first case corresponds to a rapid growth of the 
exchange coupling with decreasing interatomic distance. 
The two nearest-neighbor bonds form six local environ- 
ments with coordination numbers varying from 5 to 8 as 
shown in Fig.^b). They occur with relative frequencies 
z/^ = 17 - 12^2 ^ 2.9%, Uj^ = -41 -h 29^2 ^ 1.2%, 
uc = M- 24a/2 ^ 5.9%, u^ = -14 + 10^2 ^ 14.2%, 
j£^= 6 - 4^2 ^ 34.3%, and z^^ = -1 + ^2 ^ 41.4% 
[13 • Taking into account the short diagonals of the 
rhombic tiles increases the average coordination num- 
ber of the tiling from 4 (the value without diagonals) 
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First we discuss the Ising system. The square tile of 
the octagonal structure is non-frustrated as every pair of 
the moments can be chosen to be antiparallel (Fig.^^a)). 



If we had not taken the short diagonals of the rhombic 
tiles into account, the same would have been true for the 
entire tiling, and there would be no frustration, because 
the rhombic tiling is bipartite. Now, we consider spins 
on short diagonals as nearest neighbors, the rhombic tiles 
are always frustrated. If the energy of one nearest neigh- 
bor pair is minimized by having antiparallel spins, the 
third and forth spins cannot be chosen to minimize the 
energy of both of its neighbors (Fig. ^a)). The mag- 
netic moment will necessarily be parallel to one of the 
neighbors. For Jd < 2J two out of six possible config- 
urations have smaller energy as they possess only one 
pair of parallel nearest neighbors per rhombus instead 
of two (Fig. n^a)). In this case spins can have one of 
six possible energy values corresponding to different local 
environments (Fig. Qlb)). For Jd > 2J the four config- 
urations with two parallel bonds have lowest energy as 
their weight is smaller than that of the strong diagonal 
coupling. The second case comprises much more different 
possibilities of energy distribution. To give a quantita- 
tive description of the local frustration we introduce a 
local parameter / 
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where Ei is an actual energy of a spin i and Eid is a 
ground state energy of a relevant unfrustrated vertex. 
With this nomenclature, only the central spins of the 
vertices F and E are magnetically frustrated fp — 0.4 
and JE = 0.8 for Jd = J < 2J. The Monte-Carlo sim- 
ulations confirm our reasoning based on the analysis of 
frustration. Fig. \^ gives the frequency distribution of 
the exchange energy per atom E for two cases and a top- 
view of a portion of Ising configuration with Jd > 2 J. 
The energy distribution for Jd < 2J simply reproduces 
the frequency of 6 vertex configurations. The "up" and 
"down" configurations are perfectly ordered and coin- 
cide with the black-and-white model of Niizeki [I^. For 
large Jd we find 8 possible energy values. The "up" and 
"down" subtilings, however, are spatially disordered (see 
inset Fig.[2ti). We have calculated the magnetic structure 
factor 
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using the Monte-Carlo data for different samples. Here 
k is the wave vector and S^ is a vertical component of a 
magnetic moment at the position r. The diffraction pat- 
tern of the Niizeki configuration coincides with that of 
quantum Monte-Carlo calculations (Fig. 5c, d of Ref. p) 
and theoretical prediction 9] , while the intensity map of 
the configuration Fig.[2ti is almost structureless. It means 
that Ising solution with Jd < 2J reproduces in essence 
the antiferromagnetic superstructure, corresponding to a 
modulation vector q = (^,^,^,^)a* P ^^ the octago- 
nal tiling, whereas stronger coupling leads to a spin-glass 
state. 
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FIG. 3: (color on-line) Spatial distribution of magnetic mo- 
ments belonging to eight subtilings of a noncollinear con- 
figuration on an octagonal tiling consisting of 2193 spins. 
Jd > 2 J. The light and dark circles represent positive and 
negative x-components of the magnetization. The in-plane 
components are not given for the sake of simplicity. Average 
values of the exchange energy E and the local frustration / 
per spin are indicated. 



An exciting question is if the further minimization 
of the total energy and frustration by means of the 
noncollinear alignment of magnetic moments is possi- 
ble. At first glance the magnetic structure of the low- 
temperature pure antiferromagnetic configuration seems 
to be rather disordered. The analysis of the local en- 
ergies, however, reveals several characteristic energetic 
maxima in the frequency distribution shown in Fig.j^Ib). 
The simple existence of the peaks means that there exist 
different sorts of magnetic moments having well-defined 
relative orientation to their nearest neighbors. This ori- 
entation, however, is not associated with any absolute 
direction in space. Therefore, in accordance with the 
Mermin- Wagner theorem 14], no long-range order exists 
in two-dimensions with continuous symmetry, because 
thermal fluctuations result in a mean-square deviation 
of the spins from their equilibrium positions which in- 
creases logarithmically with the size of the system. The 
addition of a very weak anisotropy, which often exists in 
real samples, does not change the distribution of the ex- 
change energy, but permits to anchor the absolute spa- 
tial orientation of the magnetization. Nevertheless, at 
first glance the total structure still looks spin-glass like. 
In the following, we will show that the antiferromagnetic 
structure of the octagonal tiling is perfectly ordered, but 
the order is non-trivial and unusual for periodic crystals. 
We concentrate further description on 3D vector spins 
while similar results for x?/-spins have been obtained. 

To obtain an absolute symmetry axis, we apply a very 
weak out-of-plane anisotropy Ki ^ 10~^ J to the system. 
The squared vertical component of magnetization (6'^)^ 
becomes finite. The positions of the energy peaks on the 
frequency diagram remain unchanged. All maxima are 
different from those of the Ising model. It means that 
the angles between the neighboring magnetic moments 
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FIG. 4: Perspective view of a portion of a Monte- Carlo con- 
figuration on an octagonal tiling. The top view of the patch 
and the energy map are shown as insets. Magnetic moments 
are represented as cones. The cones are colored according to 
their vertical magnetization, changing gradually from red for 
"up" to blue for "down" spins. In the energy map inset, the 
colors encode the energy per moment. 



are not always equal to 180° or 0°, i.e., the magnetic 
structure is noncollinear. The different number of peaks 
— eight for Jd < 2 J and two for Jd > 2 J (Fig.I^Jb)) — 
already tells us that, in contrast to the Ising case, the 
maxima do not coincide with the 6 vertices of the tiling. 
The minimal possible local energy increases from —8 J to 
approximately —6J for J^ = J or — 5.44J for Jd = 2.2J. 
The average energy per spin, however, decreases by more 
than 0.3 J and reaches the value of ^ ^ —2. 85 J and 
E ^ — 3.30J respectively. Hence, the increase of the 
entropy permits to minimize the average frustration and 
the total energy of the system. 

Spatial arrangements of the magnetic moments as a 
function of the exchange energy are given in Fig. |31 for 
Jd < 2 J and in the inset to Fig.jSfb) for Jd > 2 J. Each 
configuration of Fig. [31 represents a certain energy range 
corresponding to one of the eight peaks in the spectrum 
of Fig. E^b). Colors represent the x-projection of the 
magnetization. The magnetic moments form 8 subtil- 
ings of different energy (£^i, . . . ^Eg) which generally do 
not coincide with a specific vertex type. The splitting 
of the energy and frustration levels is described in de- 
tail in Fig. 13 For example the vertices B and C (Fig. 
Q belong to the same energy maxima E2 but have dif- 
ferent local frustration fs = 0.24, fc = 0.11 (Fig. |S|). 
At the same time the central spin of the vertex D can 
have either the energy E^ or E4 and, therefore, can have 
two different values of the frustration f^i = 0.01 and 
fD2 = 0.11 depending on local surroundings. Thus, ev- 
ery configuration of the Fig.|31can enclose either a part of 
the atomic places belonging to one vertex type or two dif- 
ferent vertex types together. Nevertheless all structures 
have a perfect general spatial ordering. Each subtiling 
can be separated into the energetically degenerate 'right' 
and 'left' parts which also have a perfect quasiperiodic 
arrangement. However, not all 'right' or 'left' moments 
have identical orientation in space. Fig. 21 shows a per- 
spective view of a portion of typical Monte-Carlo con- 
figuration and corresponding energy map. The central 
magnetic moment has the lowest energy and belongs to 
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FIG. 5: The calculated Bragg scattering of 5*^, S^ and 5*^ 
component of magnetization for the antiferromagnetic super- 
structure. Reflexes indicated by arrows are new in comparison 
to previous studies |3]. 



the El subtiling. Its 8 nearest neighbors have identical 
energies and correspond to the energy Ej despite having 
different sets of mutual angles. The moments forming 
the next ring have energy Eq. The last ring consists of 
the alternating £^3 and Eq spins. Fig. El shows one of 
the radially symmetric vertices. However, in the octago- 
nal tiling vertices with different surrounding can also be 
found. The energy distribution is then different. Hence, 
the magnetic structure for J^ < 2 J is noncollinear and 
consists of eight interpenetrating subtilings. For J^^ > 2J 
we find only two subtilings of different energy. 

A frequency distribution of the angle between near- 
est neighboring moments shows five characteristic angles 
close to 60°, 80°, 120°, 140° and 180° for small Jd and a 
single mutual angle of 110° for large Jd. Due to this non- 
collinearity the energy of the system is decreased. The 
diffraction pattern of the whole structure is more com- 
lex than that of the Ising or the quantum-mechanical 
model. As the spin structure is noncollinear, not only 
the structure factor 5^^, but also S^^ and 5*^^ can be rec- 
ognized (see Fig.|5|. The eightfold 5*^^ and S^^ patterns 
contain additional long wave-vector peaks which could 
not be identified in the previous investigations |^. In 
dependence on the anisotropy (or on the initial random 



configuration for Ki = 0) new peaks also occur in S^^ . 
The Bragg reflexes found in our study select a subset of 
the wave vectors given in Ref. |9| where rii -\- n2 -\- ns -\- n^ 
is odd. Peaks with ni -\- 712 -\- ns -\- n^ even are extinct. 
According to the nomenclature of Ref. 9] , the following 
wave vectors can be identified: (1,0,0,0), (1,-1,1,0), 
(3,-2,1,1), (3,-1,-1,2), (1,1,-1,0), (1,0,1,-1), 
(0,2,-1,0), (0,0,1,-2), (-1,0,1,-3), (0,2,-2,1), 
(0, 1, —2, 2). Hence, the noncollinearity of the spin struc- 
ture gives rise to selection rules different from those of 
collinear models [3, [3 • With increasing sample size the 
peaks become more diffuse and may correspond to the 
diffuse scattering signal of Ref. 6] . 

In conclusion, we demonstrate that the frustrated clas- 
sical Ising system with antiferromagnetic coupling on a 
quasiperiodic octagonal tiling is perfectly ordered. All 
spins can be divided into 6 quasiperiodic (in the 3D physi- 
cal space) or 6 periodic (in 6D periodic crystal) subtilings 
of different energy. Each subtiling corresponds to the one 
of 6 vertex types of the Ammann-Beenker structure and 
is degenerated for 'up' and 'down' magnetic moments. 
Quantitatively, only two out of six subtilings are frus- 
trated with the local coefficients /e = 0.4 and fp = 0.8. 
The vector spin system admits a three-dimensional non- 
collinear magnetic structure. For Jd < 2 J, the whole 
structure can be decomposed into 8 subtilings of differ- 
ent energy which generally do not coincide with a specific 
vertex type. All subtilings are frustrated. However, the 
total degree of frustration and the energy of the system is 
minimized compared to the noncollinear case. The sub- 
tilings are degenerated with respect to the spin direction. 
The codirectional spins of every subtiling reveal perfect 
quasiperiodic ordering with a wave vector which is spe- 
cific for a given subtiling. 
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